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Rapid internal contraction boosts DNA friction 
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Macroscopic objects are usually manipulated by force and observed with light. On the 
nanoscale, however, this is often done oppositely: individual macromolecules are manipulated 
by light and monitored with force. This procedure, which is the basis of single-molecule force 
spectroscopy, has led to much of our quantitative understanding of how DNA works, and is 
now routinely applied to explore molecular structure and interactions, DNA-protein reactions 
and protein folding. Here we develop the technique further by introducing a dynamic force 
spectroscopy set-up for a non-invasive inspection of the tension dynamics in a taut strand of 
DNA. The internal contraction after a sudden release of the molecule is shown to give rise to 
a drastically enhanced viscous friction, as revealed by the slow relaxation of an attached 
colloidal tracer. Our systematic theory explains the data quantitatively and provides a 
powerful tool for the rational design of new dynamic force spectroscopy assays. 
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Life literally is a tour de force. And much of the mechanical 
load exerted on tissues and cells is supported and 
transmitted by biomolecules that can be classified as 
semiflexible polymers^ Already at physiological force levels, 
semiflexible biopolymers exhibit non-Hookean stretching 
behaviour, which sets them apart from other nanoscopic force- 
bearing elements, such as AFM tips or optical tweezers. Their 
characteristic force-extension relation under static tension^"^ has 
become the workhorse of a whole industry of single-molecule 
force spectroscopy methods, employing the polymers as gauges, 
linkers^, molecular handles^, or as a virtual magnif)^ing glass to 
explore the molecular world at ever higher spatiar and 
temporal^ resolution. However, the mechanical nonlinearity 
substantially complicates the dynamic response to external 
perturbations. If you excite an overdamped linear spring, the 
excitation always decays exponentially and the spring attains its 
equilibrium state within a characteristic relaxation time. In 
contrast, polymeric response curves exhibit no well-defined 
relaxation timescale and their relaxation generally passes 
through a multitude of power-law regimes, subtly dependent on 
the details of the applied force protocol^^'^^ Mathematically 
speaking, the response is not governed by the relaxation of a 
single eigenmode (a single equivalent spring- dashpot element) 
but by a whole hierarchy of bending fluctuations of different 
wavelengths and their corresponding relaxation times. As 
demonstrated below, such complex mode superpositions can 
give rise to somewhat unintuitive and sometimes surprising 
dynamical effects. 

A major source of complication in dynamics, compared with 
the better understood stationary situation, lies within the intricate 
coupling of the conformational dynamics of a stretched polymer 
to the solvent hydrodynamics. In its equilibrium conformation, as 
a coil of radius R, the polymer has a friction coefficient 
comparable to that of a colloidal bead of about the same radius. 
In contrast, if the polymer has a straight, rod-like conformation, it 
feels a friction force proportional to its contour length L»R. 
For a X-phdigQ DNA, LI 20, which underscores how important 
the effect can become. The result of a more quantitative 
theoretical analysis for our specific experimental set-up is 
illustrated in Fig. 1. It shows that the friction coefficient of a 
DNA-coated colloidal bead (Fig. lb) would increase by ^ 30% if 
one of the molecules of its coating layer could be permanently 
transformed into a straight conformation of length L—\6 jim 
(Fig. Ic). The purpose of the following discussion is to 
demonstrate both experimentally and theoretically that a much 
more dramatic, roughly fourfold, friction enhancement can be 
achieved if the initially stretched polymer is additionally allowed 
to recoil dynamically (Fig. Id). At the origin of this observation 
lies an excess viscous friction that is self- generated by the internal 
conformational dynamics. To our knowledge, this is the first time 
this effect has been isolated experimentally, as earlier measure- 
ments were either limited to small fluctuations around the 
stretched state^^ or suppressed the generation of excess friction 
through the attachment of large beads to the recoiling polymer 
end^^~^^. In the following, we elucidate the underlying physics of 
the enhanced friction and describe our dedicated set-up to 
measure it. We then present our measurement results and 
interpret them in the light of a systematic theory of 
nonequilibrium polymer dynamics 

In semiflexible polymer dynamics, hydrodynamic friction 
arises both due to the polymer's transverse bending undulations 
and due to the contour's longitudinal motion with respect to the 
background solvent. For a weakly bending (stretched) polymer, 
longitudinal fluctuations are geometrically strongly suppressed 
relative to transverse fluctuations. The Hnear dynamics of 
transverse bending modes therefore not only dominates moderate 




Pure DNA DNA-coated bead Hypothetical rigid-rod Our actual setup: 

scenario bead + retracting DNA 




Figure 1 | Relaxation times in a single-molecule force spectroscopy 
set-up sensitively depend on the molecular conformation. (Illustration of 
the theoretical predictions for our experimental set-up with DNA length 
/. = 16|im and persistence length ^p = 50nm, and trap stiffness 
^trap = 0.023 pN nm~^). (a) Equilibrium random coil: hydrodynamic radius 
Th ~ 0.5 |im, friction coefficient y v 6nr]r\-\ = 8.4 pN ms |j,m ~ \ and relaxation 
time T = y//ctrap~ 0.36 ms. (b) DNA-coated bead: ybead~29 pN ms |im ~\ 
T?^1.3ms (both values experimentally determined), (c) Rigid-rod 
scenario (DNA-coated bead with one DNA molecule straightened out): 
y ^ 271^1/ log(/./d) +7bead ^ 39pNms|im~\ T^1.7ms. (d) Stretch 
relaxation (DNA-coated bead with one initially stretched polymer that is 
freely retracting): 7^39 pN ms |im ~\ i^Sms (determined via 
ftrapC^^) = ftrap(0)/e). (e) The recoiling of a polymer from a stretched initial 
state is entirely driven by entropic forces or 'thermal noise' (red arrows). 
The additional contour length required to create the wiggles in the contour 
has to be pulled in from the polymer ends. The concomitant solvent flow 
(blue) exerts counteractive longitudinal drag forces that limit the maximum 
retraction speed and ultimately lead to the excess friction seen in 
scenario (d). 



transverse contour excursions, but also slow longitudinal 
fluctuations— for example, the dynamic longitudinal mean-square 
displacement of a monomer at late times t^co. But theory 
predicts that the short-time longitudinal dynamics, which 
becomes accessible by our dedicated experimental set-up, is 
strongly nonlinear and dominated by longitudinal drag. An 
example is provided by the initial stretching dynamics in response 
to a suddenly applied pulling force 

But such a fast stretching is limited by longitudinal drag for a 
relatively short initial phase of inhomogeneous tension propaga- 
tion only. Sizable longitudinal motion also occurs during the 
reverse process of contraction, when the pulled- out wrinkles 
autonomously recover to their equilibrium level after the external 
tension is suddenly released— that is, during the initial stage of 
'coiling up'^^. This 'release' scenario, sketched in Fig. le, arguably 



2 



NATURE COMMUNICATIONS | 4:1780 | DO!: 10.1038/ncomms2790 | www.nature.com/naturecommunications 
© 2013 Macmillan Publishers Limited. All rights reserved. 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms2790 



ARTICLE 




Figure 2 | Sketch of the experimental set-up. (a) One strand of DNA is 
stretched to about 90% longitudinal extension between the nanocapillary 
and the optical trap; tension is fully equilibrated, that is, fp(s, t < 0) = fp. 
(b) After turning off the external potential at time t = 0, the polynner 
retracts and the tension decays, allowing the attached bead to slowly 
relax towards its equilibrium position under the influence of the linear trap 
force ftrapCO, the polymer backbone tension fp(L,t), and the (small) friction 
fbead(0 generated by the bead itself. 

allows for the best experimental control over the comparatively 
scarcely investigated longitudinal motion limited by longitudinal 
drag, as it gives rise to an extended phase of quasistatic 
homogeneous tension relaxation limited by longitudinal 
drag^^^^ A first glimpse of this type of motion— a power-law 
growth of the retraction or 'length deficit' with t^^^ — has already 
previously been caught^^. But with our set-up, which enables us 
to directly monitor the decay of the backbone tension, we can for 
the first time inspect the underlying physical mechanism. As we 
demonstrate next, the recoiling, although it is solely driven by 
entropy, v/z., by the thermal forces resulting from the random 
bombardement of the polymer backbone with solvent molecules, 
involves such high local retraction velocities as to create the 
massive excess friction alluded to above (Fig. Id). 




-20 -10 0 10 20 30 40 50 
Time (ms) 

Figure 3 | The time-dependent force ftrapCO on the bead (deduced from 
the position of the bead in the tweezers). Apart from the (small) friction 
force fbe'ad^(f?/ktrap)7bead/5ms ^ 0.5pN onto the bead itself, ft,ap(0 can 
be identified with the force fp(L,t) exerted by the stretched polymer onto 
the bead. Experimental data were obtained from 36 independent single- 
molecule measurements (moving average over 5 ms) at five different KCI 
concentrations c = 20, 50 mM (a, 10 traces), c^lOOmM (b, 10 traces), 
c = 500mM (c, 8 traces), c = 750ml\/l (d, 8 traces). Theory curves: ftrap(0 
including the dynamic backbone tension of the recoiling polymer and the 
bead drag (solid line), using a concentration-dependent persistence 
length of ip(20 mM) = 62 nm, ^p(50 mM) = 58 nm, ^p(100 mM) = 55 nm, 
^p(500mM) = 44nm, £p(750 mM) = 41 nm (ref. 33), and a friction 
coefficient per length C\\ =0.65mPas; for comparison, the prediction for 
the hypothetical rigid-rod scenario from Fig. 1c (dotted) is also shown. 
At low salt concentrations (a) the surface charge of the capillaries has a 
significant influence on the ionic current-^^ and the observed relaxation 
characteristics exhibit a larger variability. Alternative plots with theory 
curves for ip = 50 nm (independent of the salt concentration) are shown in 
Supplementary Fig. SI. A replication of the experiment using a different 
measurement protocol is documented in Supplementary Fig. S2. 



Results 

Measurements. Our dedicated set-up to reversibly and non- 
invasively stretch and release a strand of X-phage DNA is sket- 
ched in Fig. 2. A comprehensive description of the experimental 
procedures can be found in the Methods. Briefly, a calibrated 
optical trap embedded within a microfluidic cell serves to hold or 
relocate a single DNA-coated colloidal bead^^. One of the grafted 
DNA molecules on the bead is pulled straight by applying an 
electrophoretic force with a nearby nanocapillary. After the 
electric field has abruptly been switched off at time t = 0, a high- 
speed video camera tracks the bead displacement in the trap with 
submillisecond time resolution^"*. From the bead position, we 
infer the force /ti-ap(0 exerted by the trap onto the bead. Up to the 
relatively small bead friction (that is fully included in our theory). 



/trap(0 is the value /p(L, t) of the backbone tension profile /p(s, t) 
at the grafted end s — L. Typical time traces of the force /trap(0 
exerted by the trap onto the bead are depicted in Fig. 3. 

Theory. Our systematic theory is in good quantitative agreement 
with the data (Fig. 3, lines). Since the optically trapped bead 
moves by no more than ^ 200 nm during the whole experiment, 
the bead- attached end can be thought fixed, thus closely 
approximating the "release" scenario described in refs. 16,25. Yet, 
taking also the slight bead motion into account, the theory allows 
us to reconstruct from /trap(0 the complete spatio-temporal 
evolution of the tension profile /p(s, t), including the tension 
/p(L, t) at the grafted polymer end, (Fig. 4) and the local retraction 

3 



NATURE COMMUNICATIONS | 4:1780 | DOI: 10.1038/ncomms2790 | www.nature.com/naturecommunications 

© 2013 Macmillan Publishers Limited. All rights reserved. 



ARTICLE 



NATURE COMMUNICATIONS | DPI: 10.1038/ncomms2790 




s/L 



Figure 4 | Evolution of the backbone tension along the polymer. 

Numerical solutions of equation 4 for a set of three different initial 
conditions (corresponding to boundary layer widths of about 1, 2 and 4|im 
coloured cyan, blue (dashed) and black (dotted), respectively) at various 
times t (in ms). The tension profiles quickly converge from their different 
initial states and settle onto a unique scaling function fp(s) with a slowly 
decaying amplitude, rendering the tension at the bead-laden end 
independent of the precise initial condition. 
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Figure 5 | Retraction velocities along the polymer backbone at different 
times. Theoretically reconstructed profiles v(s, t) of the local longitudinal 
contour velocity, based on the dynamic decay of the force ftrap(0 on the 
bead in Fig. 3. The computed velocity profile quickly settles down to an 
almost linear profile after the propagation layer has swept through the 
polymer (at t?^1ms). Yet, substantial retraction velocities persist near the 
free end, and the resulting enhanced friction (Fig. 1c) considerably retards 
the decay of the backbone tension at the grafted end. Inset: qualitative 
sketch (not to scale) of the recoil process and the accompanying change in 
backbone tension (from top to bottom). 



velocity along the polymer backbone (Fig. 5). In the taut initial 
state, the tension /p(s, t < 0) = /p is constant throughout the 
polymer. Once the nanocapillary voltage has been switched off, 
the polymer end is force free,/p(0, t>0) — 0. After a fast tension 
propagating period that cannot fully be resolved with our set-up, 
the backbone tension /p(s, t) approaches a characteristic invariant 
spatial shape /p(s), with slowly decaying amplitude (Fig. 4, see 
also the inset in Fig. 5 for a schematic illustration of the corre- 
sponding conformational dynamics). Notably, the force exerted 
by the polymer onto the trapped bead is predicted to be inde- 
pendent of the precise initial condition (and in particular of the 
precise value of the initial stretching force /p, in good accord with 
our data in Fig. 3). According to Newton's third law the sys- 
tematic variations in the backbone tension are balanced by the 
longitudinal viscous friction force acting along the polymer, from 



which the growth and decay of the retraction velocity along the 
polymer can be inferred^^' as depicted in Fig. 5. 

Discussion 

The exact mechanism giving rise to the observed slow relaxation 
behaviour is intuitively best understood as a dynamical equivalent 
of the static worm-like chain force-extension relation. On a 
scaling level, we balance the effective friction force C\\LdtZ (ref. 2) 
that retards the retraction of the end-to-end distance z of the 
molecule by the entropic driving force fp ^ksTL^ /[4:ip(L — z)^] 
that follows from the well-known Marko-Siggia expression for an 
almost straight conformation. Time integration then yields the 
aforementioned t^^^ — growth of the deficit L — z and the 
asymptotic power-law decay /p(t) (Cjl^A;BT/£p)^/^t"^/^ of the 
tension^^'^^. By allowing both the tension and the contraction to 
vary along the backbone, the argument can be made rigorous to 
yield the correct partial differential equation governing the 
evolution of/p(s, t) (refs 21,25). Further, including the optical trap 
and the bead friction, we obtain a slightly more complicated set of 
equations that can be solved numerically. To ease data analysis in 
future applications, we have extracted an approximate semi- 
empirical formula from this first-principles theory (the relevant 
equations of motion and the definitions of x and a are given in 
the Methods section), 

(z/«r' + 0.6(a/zrM+« l+a 



I ~ r J-./ 1 . 

This formula was validated against the exact numerical 
solutions which it matches closely over a wide range of polymer 
lengths L, bead sizes 7bead> trap stiffnesses /ctj-ap* persistence 
lengths ^p and initial values of the tension /p . It interpolates 
between the purely exponential retraction obtained for negligible 
excess friction (as, for example, in the limit of a very short/ stiff 
polymer) and the algebraic retraction /p(t) ~ that would 
result for a rigidly fixed polymer end^^ (corresponding to the 
limit of a very stiff trap). Thereby, it provides a precise 
quantification of the friction enhancement anticipated in 
Fig. 1 and a very convenient starting point for further 
quantitative applications of our dynamic force spectroscopy set- 
up. Both dimensionless parameters x ^ ^^e unambiguously 
defined in terms of experimental conditions (cf Methods), thus 
leaving no free parameter save for the effective friction coefficient 
The friction coefficient per length C|| can be estimated 
asymptotically C|| ~27i^/ log(L/a) 0.65 mPas (for a stretched 
molecule of diameter a — 2 nm (ref. 27)), and we have used this 
value for all of our theory curves. This renders the comparison 
between data and analytical theory in Fig. 3 parameter free. 

In summary, combining experiment and theory, we have 
developed a quantitative dynamic force spectroscopy set-up. 
We demonstrated that the internal dynamics of a recoiling taut 
DNA molecule generates a large excess friction that drastically 
delays the relaxation of an attached tracer bead, as quantitatively 
predicted by a simple but precise analytical approximation to the 
systematic theory. The theory also enabled us to reconstruct the 
complete spatio-temporal evolution of both the backbone tension 
and the rebound velocity along the polymer, from records of the 
tracer relaxation. In the future, it might be interesting to use this 
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Figure 6 | Analytical approximation versus exact numerical theory. The dimension less parameters a and x denote polymer size and trap stiffness, 
respectively; see equations 6 and 7. All insets show the (reduced) time-dependent force fp(L,t)/fp exerted by the polymer on the tracer bead over 
(reduced) time t/to\ numerical solution (solid lines), equation 1 (cyan dashed), pure release approximation neglecting the bead motion (blue). The 
analytical formula, equation 1, is generally accurate, with the exception of a narrow region ;(>1, 1 <a/x<8, here delineated by the dashed white line, where 
the exact solution deviates by more than O.lfp from the approximate result within the fitting interval. 



quantitative set-up to reveal how the characteristic dynamics is 
modified by the presence of DNA-protein interactions, nicks or 
supercoils^^. 

Methods 

Experimental. The microfluidic cell consists of polydimethylsiloxane (PDMS) and 
encompasses two 100 |il reservoirs connected by a nanocapillary, whose orifice size 
defines the active sensing volume for single-molecule studies. We use quartz glass 
capillaries (Hilgenberg, Germany) with an outer diameter of 0.5 mm and a wall 
thickness of 0.064 mm (20 mM) and 0.1 mm (50-750 mM), respectively. At the 
bottom the microfluidic cell is sealed with a glass cover slide of 100|j.m (20 mM)/ 
130— 160 |im (50 — 750 mM) thickness. The electric potential is applied with the 
help of two Ag/AgCl electrodes, one of them residing within the nanocapillary, 
the other located in front of the capillary tip. We use a potassium chloride (KCl) 
solution at concentrations of 20, 50, 100, 500 and 750 mM in 10 mM Tris 
buffer at pH 8. 

Our custom-built optical tweezers set-up is assembled on an optical table and 
uses a 5-W ytterbium fibre laser (YLM-5-LP, IPG Laser, Germany) at a wavelength 
of 1= 1,064 nm (ref 24). Whereas the optical trap itself is static, the microfluidic 
cell is mounted onto a xyz piezo nanopositioning system (P-517.3 and E-710.3, 
Physik Instrumente, Germany). With a range of 100 |j.m in xy- and 20 |im in 
z-direction it allows for nm-precise manoeuvring of entrapped beads in relation to 
the surrounding cell. Illumination of the region of interest inside the sample cell is 
either done using a standard white light source (DC-950 Fiber-Lite, Edmund 
Optics, USA), or an optical fibre (lOOW mercury arc lamp, LOT-Oriel, UK and 
600 |im multimode silica fibre, NA0.39, Thorlabs, UK). In combination with a 
high-speed CMOS camera (MCI 362, Mikrotron, Germany) this modular approach 



allows for real-time tracking of optically trapped colloids at up to 10,000 fps 
(frames per second) with 2 nm spatial resolution^'^. 

A Faraday cage encloses the sample cell. The externally applied electric potential 
is held constant by a commercial amplifier (Axopatch 200B, Axon Instruments, 
USA) in voltage-clamp mode. Its headstage is mounted inside the Faraday cage and 
connected to the Ag/AgCl electrodes of the microfluidic cell. Concurrent 
measurements of the ionic current through the nanocapillary^^ allow the 
simultaneous trapping of multiple DNA strands to be ruled out. 

Our DNA specimens are extracted from bacteriophage lambda (A.-DNA, 
New England Biolabs, United Kingdom). We attach biotinylated }^-DNA to 2.1 \xm 
streptavidin functionalized polystyrene colloids (Kisker, Germany)^^. With a 
binding constant of = 4 x 10~ mol the interaction between biotin and 
streptavidin is almost covalent and hence suitable for long-term experiments^ ^ 

Before each experiment we determine the nanocapillary tip diameter via its 
current-voltage characteristics^^. Afterwards, A,-DNA-coated colloids are flushed 
into the microfluidic cell and the optical trap is calibrated by analysing its 
power spectral density^'^'-^^. The proper grafting of A.-DNA to our polystyrene 
colloids is verified by repeating the power spectrum calibration for a number 
of particles. 

Tension dynamics. Apart from a short-lived initial regime (t~ns) that we cannot 
measure here, the relaxation behaviour of a freely contracting semiflexible polymer 
is governed by a generalization of the Marko-Siggia force-extension relation^ to 
space- and time-dependent backbone tension /p(s, t) 



kuTdtMs, t) 



,3/2- 



(2) 
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For a freely relaxing chain, the above equation has been derived and solved 
systematically in the literature^^'^^. Within this theoretical framework, our 
stretching apparatus constitutes an initial condition /p(s, t = 0), where the force 
gradient d/pis, t = 0) is linearly proportional to the electric field and thus to the 
resistivity determined by the nanocapillary cross section. 



and we interpolate between both extremes using the superposition ansatz with a 
mixing parameter ^, 



(1+9^2) 



1/3 



-t/T 



a/p(s,t = 0) (xE{s) oc- 



(3) 



The fairly complex distribution of forces acting within the nanocapillary 
provides a good experimental approximation to a point force applied to the 
polymer end, as details of the initially applied tension profile /p (s) quickly diffuse 
away during tension propagation, see Fig. 4. Although in principle equation 2 
breaks down close to the force-free end, and hydrodynamic interactions with 
the pore entrance or hydrodynamic boundary effects might locally increase the 
longitudinal friction coefficient these finer points turn out to be negligible 
(Supplementary Notes 1 and 2, Supplementary Figs S3 and S4). 

We now extend equation 2 to dynamic force spectroscopy assays by coupling 
one polymer end to a linear trap of stiffness Aibead and drag coefficient ybead- 
Given a certain time-dependent bead velocity Vbead(0 that must agree with the local 
polymer velocity dsfp{s = L, t)/C\\ at the bead-laden end, the dynamic tension 
profile is then determined as follows. 



/p(5,0)=/p°(5) 

/p(o,0 = o 

dsML, t)/C\\ =Vbead 



(4) 



IkBT dtfpjs, t) 



3/2 ■ 



The bead velocity Vbead is not an external parameter but depends on the 
polymeric force fp exerted on the bead. As bead and trap together constitute a 
linear overdamped subsystem with characteristic relaxation time A;trap/)'bead> 
we can infer the bead velocity from the time- dependent polymer backbone tension 
evaluated at its bead-attached end/p(L, t). 



Vbead (0 = 



-A 

7bead J 0 



(LX)e- 



Ht-i 



(5) 



This closes equation 4, which we solve self- consistently by iterating in 
Vbead(0 (Supplementary Methods). 

To obtain a uniformly valid practical approximation to the resulting force 
relaxation /p(t): =/p(L, f), we first reexpress the above system of equations in 
dimensionless form by measuring /p in units of/p, distance s in units of 
polymer length L and time t in units of 



•o^3/2 



^JkBT/ip. 



The equations of motion then read 

Jo 



where 



denotes the ratio between 'static' polymer drag and bead drag and 



^ ^11 ^trap 

^ ~ A 



4 ybead V 



(6) 



(7) 



is a measure of trap stiffness. Although the above equations may be solved 
to high accuracy using very little numerical effort (see Supplementary Methods), 
one should not underestimate the advantages of having a ready-made analytical 
expression for practical experimental work at hand. For this reason, we have 
devised the semi-empirical formula equation 1 that closely matches the result 
of the numerical integration, even in the nonasymptotic regime of initial 
decay, characterized by/p >f?{L, t) >0.1 /p. We now sketch its 
derivation. 

For /ctrap oc X ^ 00 , the bead- attached polymer end would not move at all, thus 
rendering the polymer dynamically equivalent to one half of a freely retracting 
polymer of length 2L. It has been shown before^^ that the latter scenario gives rise 
to a self-similar tension profile decaying like t~^^^ at large times; as we infer from 
our numerical data, a regularized variant /p(L, f) ~ (1 + 9f) ~ provides a 
reasonable nonasymptotic fit at the cost of a mismatching prefactor at oo . The 
opposite limit can be realized by 7bead^ oo or L^O. In both cases, the 
relaxation is completely dominated by bead friction, thus yielding an exponentially 
relaxing force /p(L, t) ~ exp( — t/T) with characteristic relaxation time 
(C||L + }'bead)/^trap which, in units of to, reads f = (1 + a)/%. For general values of a 



Next, we vary both x and a on a logarithmic scale, solving for fp numerically and 
fitting the above ansatz to the obtained solutions on the interval between /p and 

/min ~ 0. 1 /p . 

To make sure we have covered all regions of interest, we verify that our data 
closely approaches the asymptotic solutions f^ix ^ oo ) = 1 (since x^ ^ can 
be realized by taking A;trap to infinity, leading back to the pure release scenario), 
j5(a^0) = 1 (directly follows from the dimensionless equations of motion) and 
f^ix^^) = 1/(1 + oi). The latter asymptote corresponds to a virtually infinite 
persistence length » L, such that the polymer retracts as a solid object of drag 
coefficient C||L and producing purely exponential relaxation on a timescale t as 
discussed above. In contrast, the timescale of internal contraction remains 
unchanged in units of to such that tension propagation amounts to a sudden 
drop from the initially applied force/p(f = 0)//p = 1 to some smaller force 
/p(0"^ )//p = fi, which we identify with the force needed to drag along the 
polymer at an instantaneous speed i;(0+ ). As 



x{t)=Xoe 



-t/z. 



-tlx 



we have 



1 + a" 



For large a, we do not have an analytic expression for finite X- However, as for 
any finite x> ^ must approach zero as a ^ oo , we know that all nontrivial values of 
^ are realized for very large values of ^ as a diverges. The boundary condition at the 
bead-attached end thus simplifies to 



a/p(i,o 



/p(M). 



We thus only need to make sure that our numerical data approaches a self- 
similar shape /p(x, a) =/p(a/x) in the limit of large a. Taking this into account, we 
arrive at the approximate expressions for ^ and t shown in equation 1. Figure 6 
illustrates the quality of our empirical interpolation formula for nonasymptotic 
values of x and a. 
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